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Abstract 

We study the stochastic resonance (SR) of a flexible polymer sur- 
mounting a bistable-potential barrier. Due to the flexibility that can 
enhance crossing rate and change chain conformations at the barrier, 
the SR behaviors manifest many features of an entropic SR of a new 
kind, such as the power amplification peaks at optimal chain lengths 
and elastic constants as well as the optimal noise strengths. The pro- 
nounced peaks that emerge depending on the chain lengths and con- 
formation states suggest novel means of manipulating biopolymers, 
such as efficient separation methods, within undulating channels. 

The stochastic resonance (SR) is a counter-intuitive phenomenon, where 
background noise can be instrumental in enhancing coherence and resonance 
of a nonlinear system to a small periodic signal [HE]- It occurs when noise- 
induced hopping events synchronize with the signal, which itself is not suf- 
ficient to drive the system to cross over the barrier. Since the pioneering 
work of Benz et al. on recurrence of ice ages [3], the idea of SR has been 
extended to a vast range of phenomena encompassing signal and informa- 
tion processing, medicine, and biological systems [4-6]. Recently a number 
of workers reported the new mechanism of entropic stochastic resonance for 
the Brownian particle which is confined in a narrow space and thus subject 
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to free energy rather than purely energetic potential [7]. For nano-scaled soft 
matter /biological systems the ubiquitous entropic effects due to confinement 
can play an important role in noise-induced resonant effects as recently have 
been shown in the studies of SR in ion channels [H [9]. 

As interconnected, flexible systems they are, biopolymers manifest inter- 
esting cooperative dynamics in certain confined environments, external fields 
and noises. Its cooperative dynamics is important, not only in understand- 
ing how such a system self-organizes its flexible degrees of freedom, but also 
in a multitude of single molecule biophysics applications such DNA separa- 
tion, and biopolymer sequencing, etc. Nature can utilize the ambient noises 
(fluctuations) of various kinds in such biological soft-condensed matter to fa- 
cilitate the barrier crossing seemingly difficult to surmount, typically assisted 
by conformational changes. 

The study of dynamics of a polymer surmounting a potential barrier along 
with the associated SR provides a basic paradigm in which to understand the 
self-organization and cooperativity induced by the chain flexibility and fluc- 
tuations. The potential force on the chain can be traded with an entropic 
force caused by confining geometry as noted by [ID] . This make it possible 
to study the dynamics under an external potential by probing the equivalent 
dynamics within a channel with the cross section modulating over mesoscopic 
or macroscopic length scale that can be fabricated [11] . For the characteristic 
dimension of the double well potential much longer than the chain contour 
length, the crossing (Kramers) rate for a flexible chain was found to be al- 
ways higher than expected for the rigid globule it becomes in the limit of 
infinitely high elastic constant [T2l [T3] . Also the rate is found to depend on 
the conformation, coiled or stretched, which the flexible chain takes at the 
energy barrier. These behaviors suggest the emergence of an entropic SR of 
a new kind originating from chain's intrinsic flexibility, which we study in 
this Letter. 

Based on the crossing rates of such chain given in [JL 3 j , whose validity 
is confirmed by simulations [13], we recast the many-body dynamics of seg- 
ments to that of the center of mass (cm) of the chain under an effective 
potential or free energy. To identify the most salient features of this SR, 
we analyze the power amplification factor via a linear response theory. The 
power amplification factors manifest peaks as the typical signature of SR at 
optimal noise strengths, but also, remarkably, at optimal chain lengths that 
are different depending on the chain flexibility and the conformations at the 
barrier or a narrow constriction in a channel. Therefore the SR suggests novel 



2 



possibilities of understanding the biopolymers living processes, and manipu- 
lating single biopolymers utilizing its flexibility and sorting them depending 
on their lengths [15] . Earlier studies on SR of linearly coupled chains have 
mostly focused on the other extreme, i.e., the chains much longer than the 
width of the barrier. Here the barrier crossing is facilitated by a mode of 
excitation called "kink-anti kink pair", which yields the crossing rate [TB] 
and SR [T71 [TS] much different from discussed above. Because the kink is a 
localized object in the long chain limit, the activation energy is independent 
of the chain lengths and conformations, unlike in many practical cases, with 
which we are concerned. 

We consider that a linear, harmonic chain of N beads (monomers) un- 
dergoes a Brownian motion in three dimension. Assuming each bead has a 
friction coefficient 7, the dynamics of the N beads ( n— 1,2,3 . . . N) under a 
one-dimensional external potential U (x) and thermal activation is governed 
by the Langevin equation, 

7~5T = ~H 2x n - %n-l - Xn+l) - 9U Xn + £ (f) (1) 

at ox n 

where the k is the spring (elastic constant) of the chain, is assumed to 
be Gaussian and white noise satisfying 

(£n(0) = 0, (Ut)Ut + t)) = 2£> 7 5(r) (2) 

with D = ksT is the strength of the thermal noise. The external potential 
energy each bead experiences is given by 

U(*) = -&W+\0)x* (3) 

as shown in Fig. 1. The positions x = ±x m represent the two minima of the 
potential separated by a barrier of height Ub = x^Wg / 4 which is centered at 
x = 0. The parameter u B = \uJq and Uq denote the potential curvatures of 
the barrier top and the well minima respectively. 

We ask a fundamental question: How does the linearly coupled flexible 
system like a polymer chain in the double well cooperatively respond to 
a small time-periodic force in the presence of the ambient thermal noise? 
The determinant parameters are the noise strength (temperature), the chain 
length, and the elastic constant, given the double well potential parameters 
as fixed. By such coupling the SR is enhanced above that of a single bead 
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U(x) 




Figure 1: Schematic diagram for a flexible polymer chain in a symmetric 
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double-well potential U(x) = — C^-)x 2 + j(^-)x 4: . The potential wells and 
the barrier top are located at x = ±x m and at x = 0. We consider the case 
of the potential where NUb 3> -D, and L ^> X m so that the potential is 
designed to vary on a mesoscopic or macroscopic scale. 



as already has been studied in earlier studies [17,18]. In the globular limit 
where the coupling (spring) constant is infinity, the polymer moves as a 
single Brownian particle (a rigid globule) in the bistable potential, whose SR 
behavior is well-known. Now turning on the internal degrees of freedom that 
give rise to the chain conformational flexibility and variability, how will the 
SR be affected? 

We base our theoretical development of SR on previous theoretical results 
[13] for the rate of the chain crossing over the potential barrier by thermal 
activation in the absence of driving force, which is found to be consistent 
with numerical simulations [13, 14]. The rate is given by 



AF 



R = R e~— (4) 
where 

Here the Rq is the rate in the globular limit, k — > oo. As the chain attains 

AF 

the flexibility by decreasing k, the rate is enhanced by the factor e ~ where 
AF is the free energy change of the chain during the barrier crossing, which 
depends on temperature. 
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It is because the chain, a flexible extended object as it is, experiences an 
entropy increase in surpassing the potential, i.e., from the confining well to 
the unstable barrier top, thereby reducing the activation energy barrier [13]. 
AF depends on D as shown below. Depending upon the chain length TV 
and the rescaled spring constant k = k/u B , the flexible chain takes either 
the coiled or stretched conformation at the barrier top (at the transition 
state). For the regime of parameters satisfying q > where 

the chain retains the coiled conformation, the free energy change is expressed 
by 



AF = -Din 





1/4 j 













sin(NyJu%/k) 




-TAS C (7) 



and f(a) = ^ JZo d Qe~ {a/2)Q2 ~ 3/8Q4 with a = 2g(^) 5 . It should be 

noted that AF' is given solely from the entropy change AS C mentioned above. 
On the other hand, when q < 0, the chain takes a stretched conformation at 
the barrier top, thereby reducing the energy barrier by the amount 

- W ,76 - D In \V2 ( ^ | | ^ | V4 ( g(a) 




^sin(N^oJ B 2 /k) i 

for small values of \q\. In Eq. (9) g(a) = ^e~ a2 / 6 dQe~^ Q2 - 3 / 8Qi 

where (cJs) 2 = (1 — 2\q\)u B . Eq (8) includes the decrease in the internal 
energy due to stretching ( the first term) as well as the entropy increase 
(the second term) in crossing the barrier in stretched state. The transition 
from coiled to stretched conformation at the barrier are incurred as the chain 
length increases above N c or the spring constant decreases below a critical 
value k c with the potential parameter fixed. 

Driven by a force F(t) = F cos (fit) on each bead, where F Q and fl are 
the amplitude and angular frequency, the center of mass position X of the 
chain, the reaction coordinate for the barrier crossing, evolves following the 
Langevin equation under an effective potential or free energy function F(X) 

= -^ffl + NF cos(Qt) + C(t) (9) 
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where £ (t) represents the random force characterized by (£ (£)) = 0, (£ (t + r)) 
2DNjS{t). 

It is reasonable to assume the free energy F(X) associated with the cm 
position X takes the form similar to Eq. (3), 

F(X) = -^X 2 + ||x 4 (10) 

with the barrier curvature Q B and height AF = ^Q B x"^. The Kramer rate 
for the cm regarded effectively as a single Brownian particle is 

F = — — e o (11) 

which is equated with Eq. (4) to yield 

AF 9* NU B + AF' (12) 

and thus 

i— ( 4AF' \* 

n " = i 1 + «=S*J ' (13) 

both of which depend on JV as well as the D via AF . The free energy 
function, rather than purely energetic potential, under which the particle 
representative of the chain is moving, signals an occurrence of an entropic 
stochastic resonance [7]. 

How can we realize the dynamics of the chain under the one-dimensional 
free energy F(X) in a manageable experiment? To this end we convert 
the dynamics to that of the chain flowing within the channel of undulating 
cross section A(X). Since we envision the chain as a single particle moving 
subject to F(X), such a conversion is possible if F(X) = —DlnA(X), or 
A(X) = Aoexp(—F(X)/D, where Aq is the cross section at the point where 
F(X) = 0. Therefore the crossing dynamics and the associated SR can 
be achieved for a chain within such a fluidic channel subject to an external 
periodic forcing, for example a single stranded DNA within a channel subject 
to AC field. 

To gain understanding of the major salient features of the distinctive 
nature of the polymer SR, it suffices to consider the linear response of the 
chain cm to the small driving forces. In response to the weak driving force, 
the average cm position is given by 



X(t) = NFo\x(£l)\cos(nt-(/)) (14) 
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where x(fy ^ s the Fourier transform of the response function [Tj5] given by 

X(^) = + *x"(0) = D~lMin and ^ is the phase delay given by 

tan -1 ^//m while the coherence of the system is measured by the real part 
of x(O), the SR intensity is quantified by the power amplification factor (the 
ratio of power stored in the response of the system to the power of the driving 
force with frequency): 

where (X 2 ) = JX^'^dX/ J e'^dX. 

The r\ indeed shows nonmonotonic noise strength dependence manifest- 
ing a peak at an optimal noise strength, which is the typical signature of 
SR. Via dependence of the rate R and (X 2 ) on the conformational variabil- 
ity, the SR for the flexible polymer manifests distinctive dependence on the 
noise strength and additional dependence on the chain length, characteristic 
of a polymer entropic SR. Below we quantitatively discuss the power ampli- 
fication factor 77 for the polymer conformation both below and above coil-to- 
stretch transition. We introduce dimensionless parameters: D = D/ugX^, 
Cl = '-fQ/uj 2 3 and k = k/uj 2 . From now on all the quantities are rescaled 
(dimensionless) so the bars will be dropped. 

Let us first calculate 77 as a function of rescaled noise strength D for the 
rescaled spring constant k = 200 which can be reasonable for a real flexible 
polymer if lo 2 b is chosen properly. The Figures 2a and 2b depict the rescaled 
77 as a function of noise strength D for Q = 10~ 5 and Q = 1, respectively 
for various transition states of polymer conformation; the coiled state with 
N = 30 and the stretched state with N = 50, both of which are compared to 
globular state (the case of infinite k) with N = 30. Hereafter in the figures, 
the coiled, stretched and globular states are depicted as dashed, solid and 
dotted lines. As shown in the Figs. 2a and 2b, for all cases, the 77 reach 
their maximum value 77^ at certain optimal values of D, denoted by Dr. 
When the chain is in coiled state, 77^ is much larger that of globular state. 
This signifies that its intrinsic flexibility facilitating faster crossing allows the 
polymer chain to respond in a more cooperative and coherent manner to the 
external signal. 

For very small Q, as depicted in Fig. 2a the optimal noise strength 
corresponding to the peaks satisfy D RC < D RG < D RS , where the subscripts 
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Figure 2: (a) 77 versus D for Q = 10 -5 . (b) 77 x 10 4 versus D for Q = 1. The 
solid line stands for stretched state (N = 50), the dashed line for coiled state 
( N = 30) at the barrier top and dotted line for the globule (N = 30) for 
parameter value k = 200 for which Nc = 44.4 is the critical length for the 
coil-to-stretch transition. 



C, G and S denote the three conformational states. To understand this, we 
note that the time scale matching condition between the mean crossing time 
and half of the period of oscillation, R ^ R ^ = which provides a reasonable 
approximation for the resonant condition. With R{Dr) Pd 3 x 10 -6 given 
thus, and from the R obtained from Eqs. (7) and (8) (see Fig. 2a), the 
above inequality indeed holds. The figure 2a shows that r] RC > Vrg > Vrs 
which can be understood as follows. At very low frequents, the rj, Eq. (15) is 




which implies that the resonance mechanism 



is governed by chain static susceptibility |x(0)| = - D - of the cm in response 
to a constant force. Finding the static susceptibility for the D R of each 
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transition state (D RC < D RG < D RS ), one indeed affirms the inequality for 
the r] R . 
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Figure 3: (a) Polymer crossing rate R versus (times 10 6 ) D for small values 
of D. The dotted, dashed and solid lines stand for globular {N = 30 and 
N = 50), coiled (N = 30) and stretched state (N = 50) at the top of the 
barrier, (b) Polymer crossing rate R versus D for large values of D. 

The situation is different for large driving frequency where the resonance 
conditions are met for larger D R (see Fig. 2b). While the inequality for D R 
is unchanged, that of r] R is not; r] R for stretched state is highest. This is 
understandable if we note for large Q the rj is largely governed by the rate, 
rather than the static susceptibility. As the chain length becomes longer than 
the critical value for this large D, this flexible chain can extend to lower the 
free energy barrier and enhance the crossing rate, as shown in Fig. 2b. This 
results in the stretched conformation enhancing SR (solid line) than the coiled 
or globular states at these noise strengths as the Eq.(20) indicates. This 
means that the chain capable of the conformational transition can utilize such 
flexibility to enhance the cooperatively and coherence to external influences. 
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Figure 4: (a) r] as a function of N for noise intensity D — 1 (the top curve) 
and D = 2 (lower curve), (b) 77 (times 10 2 ) as a function of iV for noise 
intensity D = 4. Other parameters are fixed as Q = 0.01 and fc/cu^ = 200. 

Another aspect of this entropic SR is the phenomenon above is shown in 
the Figures 4a and 4b which depict rj a function of chain length with k = 200 
and Q = 10 -2 for three different noise strengths, D — 1,2, and 4. For D — 1, 
the chain manifest the response of coiled conformation for all values of N 
which are mostly smaller than iV c = 44.4. For D = 2, the chain responses 
in coiled conformation for N < N c but responses in stretched conformation 
for N > N c . For D = 4 the chain manifest smaller rj for all chain of length 
N > N c which are all in stretched conformation at the barrier top. Each 
of these lines have a peak at a certain optimal 7V which is increasing with 
D. Even with the temperature fixed, as is usual in biological systems, a long 
chain by stretching can escape the barrier or the narrowest constriction in a 
most coherent and resonant way to an external signal. 
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Figure 5: (a) rj versus D. The solid line represents k = 0.8k c (stretched 
state), dashed line k = 1.5k c (coiled state), and dotted line k = oo (globular 
state), (b) 77 versus k for N = 50 and D — 1.5. In both figures we fix 

n = 0.01. 

Let us now fix the polymer length N = 50 and vary the rescaled spring 
constant k where k c = 253.3 is the critical spring constant that demarcates 
the polymer conformations at the barrier. Figure 5a exhibits 77 versus D, 
where the curves depict k = 0.6k c , k = 2.5k c and k = 00, respectively 
representative of stretched, coiled, and globular conformation at the barrier. 
As the spring constant (stiffness) increases, the chain shows less resonant 
peaks at larger optimal noise strength. This means that larger thermal energy 
is required to drive the stiffer chain to move coherently with an external 
signal. On the other hand Figure 5b exhibits the plot of rj as a function of 
k for N = 50. As k becomes very small, 77 tends to be very small as the 
whole monomers becomes non-interacting (non-cooperative). When k — > 00, 
77 becomes small again as globular polymer lacks the enough flexibility to 
conform to external driving. At certain optimal value, in between 77 attains 
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an optimum value which is an increasing function of D. 

In summary, we have studied the stochastic resonance of a single flexible 
polymer moving in a bistable potential, or flowing within an equivalent fluidic 
channel, both of which modulate over a mesoscopic or even a macroscopic 
scale. At an optimal noise strength, the dynamics of the chain, otherwise 
random, shows a SR, i.e., it moves in coherence and resonance with a periodic 
driving force. Even with the noise strength (temperature) fixed, the polymer 
display a novel kind of entropic SR by responding cooperatively to external 
driving in a most efficient way at optimal chain lengths and elastic constants, 
owing to chain flexibility and conformational transition. Utilizing their self- 
organizing behaviors, we may learn bio-molecular machineries of living, and 
clever ways of manipulating them such as efficient separation methods. 
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